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^ ! Abstract 

^ I A mono-component is a real-valued signal of finite energy that has non-negative in- 

Q ' stantaneous frequencies, which may be defined as the derivative of the phase function of 

the given real-valued signal through the approach of canonical amplitude-phase modula- 
tion. We study in this article how the amplitude is determined by its phase in a canonical 
amplitude-phase modulation. Our finding is that such an amplitude can be perfectly re- 
constructed by a sampling formula using the so-called generalized sampling functions and 
their Hilbert transforms. The regularity of such an amplitude is identified to be at least 
O I continuous. Meanwhile, we also make a very interesting and new characterization of the 

band-limited functions. 
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^ 1 Introduction 
(N 

Any real- valued non- stationary signal / of finite energy, that is, / is in the space L^(M) of square 
integrable functions on the set M of real numbers, may be represented as an amplitude-phase 
modulation with a time-varying amplitude p and a time-varying phase (j) where phase (j) is, in 
general, nonlinear. Specifically, the value of / at t G M may be represented as 

f{t) = p{t) cos <j){t). (1.1) 
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Unfortunately, this type of representation is not unique because the modulation is obtained 
through a complex signal that can have various choices of the imaginary part. However, one 
can determine a unique such factorization ( II. ip by using the approach of analytic signals. 
Indeed, let A{f) be the analytic signal associated with f with the characteristic property 

= 'J'^^" (1.2) 

10 if w < 0, 

where for any signal g G L^(M), g = J^g is the Fourier transform of g defined at ,^ G M by the 
equation 

m = {^gm--=^ [ 9it)e-'^'dt. (1.3) 
Equation (11. 2p is equivalent to for t G M, 

A{fm = f{t) + mf{t), 

where the operator H : L^(M) — L^(M) stands for the Hilbert transform, and for / G L^(M), 
"H/ at t G M is defined through the principal value integral 

nfit):=p.v.- [ /Mdx = lim [ ^dx. 



The value A{f){t) at t G M is complex which can be written into the quadrature form 

AU){t) = p{ty^^'\ 

Under the conditions that the derivative value 0'(t) is non-negative or non-positive for all 
t G M, the quantities p{t) and 0'(t) are called the instantaneous amplitude and instantaneous 
frequency at t G M, of the real signal /, respectively. The corresponding modulation (II. ip is then 
called the canonical amplitude-phase modulation, or canonical modulation for short. The signal 
/ with such defined non-negative instantaneous frequencies is thus called a mono- component. 
A large body of hterature addresses this problem, see for example, [U [21 El [HI IHl HI 116] . 

Constructing the canonical pair (p, 0) of the instantaneous amplitude and phase is important 
in the theories of analytic signals. It is equivalent to the problem of seeking the function pair 
(p, 0) such that for t G M, the following equation holds true 

?/(p(-)cos0(-))(t)=p(t)sin0(t). (1.4) 

We remark that (II. 4p can be apparently considered as a special case of the Bedrosian identity 

nfg) = f-Hig). 

In [1], the author proved that, if both f,g belong to i^^(M), / is of lower frequency, g is of 
higher frequency and f,g have no overlapping frequency, then T-i^fg) = fHi^g). This classic 
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result of Bedrosian is not useful for constructing a mono-component. The reason lies in that 
the requirement of both / and g in L^(M) is invalid. 

Recently, an important phase function that renders mono-components was given in |13j . 
The phase function is defined through the boundary values of a Blaschke product on a unit 
disk A := {z : z G C, \z\ < 1}, where C indicates the set of complex numbers. Specifically, for 
a G (—1, 1), the Blaschke product at z G C\{l/a} is given by 

= i^. (1.5) 

Subsequently the non-linear phase function, denoted by 6a, is defined at t G M by the equation 

e^'-W .= Baie"'). (1.6) 
If we recall that the periodic Poisson kernel Pa whose value at t G M is given by 



Pait) := ^ " (1.7) 

1 — 2a cost + 



then by taking the derivative of both sides of equation f ll.6p . we find that the phase 9a is an 
anti-derivative of pa, and its derivative is always positive, that is, 

^^a(t)=Pa(t)>0. 

We shall in this paper characterize the amplitude function p of finite energy when the phase 
function (j) is chosen at t G M by 



(f){t)=ea{t) = / Pa{x)dx 
J[0,t] 

such that equation (11.41) is satisfied. Our main result indicates that such kind of amplitude 
can be perfectly reconstructed in terms of a sampling formula using the generalized sampling 
function whose value at t G M is given by 

. /,x sin^„(t) 

smca(t) := . (1.8) 

In Section 2, we review the construction of the generalized sampling function and discuss 
some properties pertained to it. In Section 3, we introduce the concept of Bedrosian subspace 
of the Hilbert transform and investigate some properties of functions in this space. In Section 
4, we make an important observation when a linear phase is chosen, the amplitude function 
must be bandlimited in order to satisfy equation (II. 4p . In Section 5, we present our main result 
in Theorem 15.71 
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2 Generalized sampling functions 



Not very surprisingly the function sinCa has many properties that are similar to the classic sine 
defined at t G M by the equation 

sint 

smc()f:) := 

Those properties include cardinality, orthogonality, decaying rate, among others. In the special 
case a = 0, the function sinca reduces to the classic sine, which will become clear later. Let us 
first review the approach to obtain an explicit form of sinca. 

The classic sine function is fundamentally significant in digital signal processing due to the 
Shannon sampling theorem [T31 [151 E| • The Shannon sampling theorem enables to reconstruct 
a bandlimited signal from shifts of sine functions weighted by the uniformly spaced samples of 
that signal. Recently efforts have been made to extend the classic sine to generalized sampling 
functions, for example, in O El [7]. Intuitively, the spectrum of the sine function is just the 
indicator function of a symmetric interval of finite measure. Hence, authors in p] are inspired 
to consider functions with piecewise polynomial spectra to replace the usual sine function for 
the purpose of sampling non-bandlimited signals. One kind of generalized sampling functions 
given in [5], denoted by sinc^ that is related to a constant a G (—1, 1), is defined as the inverse 
Fourier transform of a so-called symmetric cascade filter, denoted by Ha- Specifically, 

sinc„ := + a)T-^Ha. (2.9) 

Let N be the set of natural numbers, Z be the set of integers, and Z+ := {0} U N. Let X 
be a subset of M, and for g G N, we say a function / is in L'^{X) if and only if the L'^{X) norm 



q,X 



1/9 

\f{t)\'dt] <oo. 



X 



Similarly, let Z he a subset of Z, a sequence y := {yk : G is said to be in if and only 
if the l'^{Z) norm 

\kez J 

The symmetric cascade filter Ha is a piecewise constant function whose value at ^ G M is 
given by 

Ha{0:= $^a"X/„(0, (2.10) 

where xi is the indicator function of the set /, and the interval J„, n G Z4., is the union of two 
symmetric intervals given by the equation 

J„:=(-(n + l),-n]U[n, (n + 1)). 
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Of course, we have that H G L^(M) n Z/^(M) because the sequence h := (a" : n G Z+) G 
n/^(Z_|_), and hence sinc^ G L^(R) since the Fourier operator is closed in L^(]R) and sinc^ 
is continuous because Ha G L^(M). 

The symmetric cascade filter Ha can be associated with an analytic function F on the open 
unit disk A defined at z G A by 

F{z) ■= a^'z''. (2.11) 

Thus, by substituting equation fl2.10p into equation (12. 9 p and making use of equation (12. lip an 
alternative form of sinca(t), t G M in terms of F can be found as 

sinca(t) = (1 + a)sinc Re |F(e")e^"| , (2.12) 

where Re(z) is the real part of a complex number z. 

A very interesting fact, as discovered in the paper [3], is that the function F is linked to the 
Blashke Product Ba by the equation 

(l + a)F(^):=^?^^4r^- (2.13) 

Plugging the formula (I2.13P into equation (I2.12p we readily obtain the explicit expression of 
sinca given earlier in equation (II. 8p . 

The next two formulas shall be used later. Expanding the left-hand side of equation (II. 6p 
using Euler's formula and separating the real part from the imaginary part of the right-hand 
side, one obtains that for t G M 

^ / X (1 — a^) sint , , , , 

sin^„ t = ^——^ = pa{t) sint 2.14 

1 — 2a cost + 

and 

^ , , (1 + a'^) cost — 2a . _ ^, 

^osOait) = \ ' 2.15 
1 — 2a cost + a^ 

We next list some properties of the function sinc^. 

Proposition 2.1 Let the generalized sampling function sinCa be defined by equation (I2.12p or 
equation (II. 8p . Then the following statements hold. 



1. Fort G R, 

(1 — a^) sint 
2a cost + a^ t 



= l-2acost + a^ ^ = J>a(t)sinc(t). (2.16) 



2. 



Tsmca= J^il + a)Ha. (2.17) 
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3. sinca(n7r) = jrf^n.O; where 5nfl = 1 if n = and Snfi = if n E Z \ {0}. 

4- sinCa is even, bounded, infinitely differentiahle. 

5. |smca(t)| < Y^Yq^ /or t G M, and sinc^ G L'^iR). 

6. The set {sinca(- — mi) : n G Z} is an orthogonal set, that is 

I ■ ■ ( w (1 + X 
(sinca,smca(- -nvr)) = 7r^„,o, 

1 — a 

•where {u,v) = J^u{t)v*{t) dt denotes the usual inner product of the two functions u,v E 
L^(M), and v* is the complex conjugate of v. 

Proof: Equation fl2.16p directly follows by substituting equation fl2.14p into equation (ll.Sp . 
Equation fl2.17p is obtained by taking the inverse Fourier transform of both sides of equation 
fl2.9p . The third and fourth statements directly follow from equation (I2.16p . The fifth statement 
follows from equation fl2.16p and noticing sinc(t) < and Pa{t) < jz^ for any t G M. The last 
statement is a special case of Corollary 3.2 of [5]. For the convenience of readers, we provide a 
direct proof here. By Parseval's theorem and equation fl2.17p we have 




where, in the last equality we have used the orthogonality identity Jj^ gi^'ra: _ 2Snfi, k G 
The interchange of the integral operator and the infinite sum is guaranteed by the absolute 
convergence of the series. □ 

In figure [2?n we show an example of the Fourier transform pair sinca and y^(l + a)iJa with 
a = 0.5. In the plot of sinCa, the graph of the standard sine is also shown, which corresponds 
to the case a = 0. 

Next we review several basic properties of the Hilbert transform which we will need fre- 
quently later. These properties can be found, for example, in the book [10]. First the Hilbert 
transform is an anti-involution, that is, 

n"^ = -X (2.18) 
where X is the identity operator. Second the operator Ti is anti-self adjoint, that is, 

{nu,v) = {u,-nv). (2.19) 
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Figure 2.1: Left: the graph of ^^(l + Ci)Ha, a = 0.5; Right: the graph of sinCa, a = 0.5. 

Third, for any / G L^(]R) and t G M, the composition of the Fourier transform and the Hilbert 
transform is given by 

J-(H/)(t) = -isgn(t)^/(t), (2.20) 

where sgn(-) is the signum function having values defined by sgn(x) = 1 if a; G IR+ := {t G M : 
t > 0}, sgn(x) = -1 if x G M_ := {t G M : t < 0}, and sgn(O) = 0. 

Theorem 2.2 The system 

$ := {sinca(- - 2kTr),nsmCa{- - 2kTr) : k e Z} (2.21) 

is an orthogonal system in L^(M). 

Proof: By the third statement of Proposition 12. H we have that 

1 + a 

(smca,smca(- - 2k7r)} = dok- 

1 — a 

Invoking equations fl2.18p and fl2.19p yields 

("HsinCa, 'HsinCa(- — 2kn)) = (sinca, — ?{^sinca(- — 2k7i)) 

= (sinca,sinca(- - 2k7r)) = ^-^6ok. 

By Parseval's theorem and equation (12.171) we have 

/ sinCa(t)HsinCa(t - 2A;7r) dt = + ^)'^^ j Hl{t)sgn{ty^''^' dt. (2.22) 
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Invoking the expression (12.101) of Ha, equation fl2.22p becomes 

sinCa(i(:)'HsinCa(t — 2A;7r) dt 
;i + a)27ri 



2 

[l + afm 
2 

'l + a)27ri 



J2 ci'-XiAt>gn{t)e'''-' dt 



2n 



,2n 



ne2 



sgn(t)e'''="* dt 



gi2fc^t dt 
.J [n,n+l) 



dt 



(-(n+l) -n] 



where the interchange of the integral and the infinite sum in the second equahty is guaranteed 
by the absolute convergence of the series. 

When k = 0, the difference in the pair of brackets is zero, while when A; G Z \ {0}, each 
integral inside the pair of brackets is zero. Therefore we have 



(sinca, ('Hsinca)(- — 2A;7r)) = 



for any G Z. 



□ 



3 Bedrosian subspace of the Hilbert transform 

In this section, we pay attention to the set 

5a :={/:/ G L\R), cos ^a(-)) = /sin^,(-)}. (3.23) 

It is clear that the set Sa is a subspace of L^(M) due to the linearity of the Hilbert transform. 
The subspace Sa shall be called the Bedrosian subspace of the Hilbert transform. We first make 
a simple observation that we shall need frequently later. 

Lemma 3.1 Equation (11.41) is true if and only if for t G M, 

n (p(-)e''^-) it) = -ip(t)e'^*. (3.24) 

Proof: Applying the Hilbert transform "H to both sides of equation (ll.4p and utilizing equation 
(|2:T8|) yields 

nipi-) sm{(f)-)){t) = -p{t) cos(0t). (3.25) 

Combining equations (II. 4p and (I3.25P produces equation (13.240 . Thus equation (ll.4p is equiv- 
alent to equation (13.241) . □ 

Lemmas I3.2H3.4I are several technical lemmas we need in the sequel. 
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Lemma 3.2 Suppose that f G Z/^(M) and g G Z/^(— |, |) is a t -periodic function having Fourier 
series g{t) = ^^c^e^'^"* with Ck = ^ g{t)e~^~^^ dt. Moreover, the series ^^c^ is 

absolutely convergent. Then fg G L^(]R) anc? zis Fourier transform at ^ eM. is given by 

Proof: The fact that g G L^(— |, |) imphes that \g{t)\ < c, for a.e. t G M and some constant c. 
Thus f'^{t)g^{t) dt < f^{t) dt < oo. That is, fg e L'^{R). By the definition ([TS]) of the 
Fourier transform, we have 



ifgnO = [ f{t)(Y,Cke'''^']e-''^dt 
= $:c.//(t)e-*(^-'=^)dt 



z — / J- 



where the interchange of the integral and the infinite sum is guaranteed by the absolute con- 
vergence of the series ^^c^. □ 



We say that a complex r-periodic signal g is circularly analytic if its Fourier expansion is 
the one-sided series g{t) = c^e'"'^*, t G M. 

fcez+ 

Lemma 3.3 Suppose that f G L^(R) is a real-valued function and the complex function g is 
T-periodic and circularly analytic with real Fourier coefficients Ck, k G Z_|_. Then 

U{fg){t) = -if{t)g{t) (3.26) 

if and only if for a.e. ^ G M_, 

^kh^ - -k) = 0. (3.27) 

Proof: Equation (13.261) indicates that fg is an analytic signal. This is equivalent to say that 
fg has no negative frequency, that is, {fg)^{0 = 0? ^ ^ By Lemma [3.21 and the fact that 
g is circularly analytic, we obtain equation (I3.27p . □ 

Lemma 3.4 The 2it -periodic function g := e'^"*-'-* is circularly analytic and the Fourier expan- 
sion o/e'^"*^'-' at t eM. is given by 

e^eait) = _a + ^(i _ a^)^^-!^!*^*, (3,28) 

fceN 
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Proof: Recalling equation ( II .Gp . the Fourier coefficients Ck, k E 'Z, of g is given by 



Ck = i^ 1 dt 

2vr 1 - ae'* 



Let R.ez{f{z); c) indicate the residue of the function / at the point c G C. Denote the boundary 
of the unit disc A by OA. For = 0, we have Cn = 7^ z^dz ^ ( ; ) = -a. In 

the case is a negative integer, by using the Cauchy theorem we get that 



27r ./(_7r,7r) 1 ~ '^e^* 27ri /^^ 1 — 



because the integrand is analytic in the unit disk A as |a| < 1. We are left to consider the case 
of G N. Using the residue theorem and the formula for j G N, 



d-' f z — a 



dz^ \1 — az 

leads to when k E N, 



[l-a^)a^'^j\{l-az^-^-^ 



= 2^ £ (1 - a.)" = ( (1 -az)z^^^ ' ° 



1 d''^ / z — G 

77 lim -— ) = fl - a^)a 

k\ z^o dz" \1 — az 



2\ „fc-l 



The proof of this lemma is completed. □ 

Next proposition indicates that the subspace Sa is invariant under the Hilbert transform. 
Proposition 3.5 Suppose that p G Sa, then Hp G Sa- 

Proof: We first claim that the complex- valued signal {p + i'Hp)e^^° is an analytic signal. Indeed, 
by Lemma (3.41 the signal e^^" is circularly analytic, and we write e^^"'-*^ = 'Ylik&+ 9k^^^^ ■ Then 
by Lemma 13.21 we have for G M, 

^((p + iHp)e^^")(0 = J]^7fc-^(p + iW(e-fc). (3.29) 

Note that (p + iHp) is an analytic signal thus it has no negative spectrum, that is J-" (p + iHp) {C) 
when ^ G M_. Consequently, when ^ G ]R_, 

Y.QkJ'ip + mp) (e-A;) = 0, 

fcez+ 

which implies by Lemma 13.31 that 

n [(p + i?{p)e'^"] = -i(p + mpy'^\ 
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Simplifying both sides leads to T-L[{'Hp)&^''] = — i[(?{p)e^^°], which indicates that Tip G Sa- □ 

The next lemma indicates that any real-valued function p & Sa is completely determined 
by its spectrum px{-i,i) on the interval (—1, 1). 

Lemma 3.6 The real-valued signal p ^ Sa if and only if for a.e. G M \ {0}, n G 

p{^ + sgn{On) = a-m- (3-30) 

Proof: By Lemma 13.11 p G iSa is equivalent to for t G M, 

n (p(-)e''"") (t) = -ip(t)e'^''W, (3.31) 

which is equivalent to saying that the signal p{-)e^^'^^'^ is an analytic signal. Applying Lemmas 
13.31 and [3^ to the signal p(-)e'^"^'^ yields that p E Sa holds if and only if for a.e. ^ G M-, 

ap(0 = (l-a2)5^a^-W-^)- (3-32) 

fceN 

In equation fl3.32p . separating the first term from the rest in the sum on the right-hand side, 
we have 

am = (1 - a') ( p(e - 1) + E ^'''p(^ - 

= p(e - 1) - a2p(e - 1) + (1 - a') J2 - k). (3.33) 



feel- 



On the other hand, multiplying both sides of equation fl3.32p by a and replacing ,^ by ^ — 1 we 
also have 



aV(e-l) = {l-a')J2a'p{^-l-k) 

ken 

= {1-a') J2 a'-'p{^-k). (3.34) 
fcei+N 

Combining equations (13.331) and (13.341) deduces that for a.e. ^ G M-, 

p(e-l) = ap(0. 
By induction, we can get that for n G Z_(. and a.e. ^ G M_, 

p(e-n) = a"p(0- (3.35) 
Reversing the above calculations reveals that condition (I3.35P is also sufficient for (13.320 . 
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Taking the conjugate of both sides of equation (13.320 and using the Hermitian property of 
p, we obtain an equivalent equation 

am = (1 - «') E ^'"''^(^ + ^) (3.36) 

for a.e. ^ G ]R_. A similar discussion to the case ^ G M_ gives that equation (13.361) is equivalent 
to 

p{^ + n) = a^m (3.37) 

for n G Z+ and a.e. C, G Finally combining equations (I3.35P and (I3.37P produces (13. 30 p . 
□ 



Proposition 3.7 

sinCfl G Sa and "HsinCa G Sa- 

Proof: By Lemma 15751 it suffices to show sinCa G Sa- Observing the identity Ha{-+n) = a^"^Ha{-] 
for a.e. ^ G M \ {0}, n G Z and recall equation (12.170 we thus have 

J-(sinc,)(e+|n|sgn(0) = al"l^(sinc„)(0. 
Hence by Lemma [3761 we conclude that sinca G Sa- 



□ 



4 An observation— how a linear phase determines the am- 
phtude 

In this section, we specifically consider the case when the phase in equation (II. 4p is a linear 
phase and investigate the representation of the corresponding amplitude. The result is given in 
the following theorem. We use the notation supp(/) for the set of real numbers on which the 
values f{x) at x G M are nonzero. 

Theorem 4.1 Suppose that 'j is a positive real number and p is a non-zero real signal in L^(M) . 
Then the following equation 

?^(p(-)cos(7-))W=pWsin(7t) (4.1) 
holds if and only if p is handlimited with supp(p) C [—7,7]. 
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Proof: Lemma [3.11 implies that equation (14. ip is equivalent to 



n (p(-)e'^-) (t) = -ip(t)e'^*, (4.2) 

or equivalently, the complex signal p{t)e^"'^ is an analytic signal. Note that the Fourier transform 
of the function on the left-hand side of equation (14. 2 p is given by — i sgn(^)p(^ — 7) and the 
Fourier transform of the right-hand side of equation (14.21) is —'ip{^ — 7). This leads to the 
equivalent equation of (14. ip in the frequency domain 

(sgn(O-l)p(e-7) = 0, (4.3) 

i.e. 

p(O=0, for ^G(-oo,-7]. (4.4) 

Thus if p is bandlimited with suppp C [7, 7], then clearly equation (14. 4p is true, hence equation 
(14. ip is true. On the other hand, assuming that equation (14. ip is true, we obtain that p(^) = 
for ^ G (—00, —7]. Since p is real- valued, by the Hermitian property of the Fourier transform 
of a real signal, p(^) = for ^ G [7, -t-cx)). Consequently we conclude that p is a bandlimited 
function with its support in the frequency domain belongs to [—7,7]. □ 

Remark: From this theorem we know that the amplitude function p can be represented by 
shifts of sine function if and only if equation (14. ip is true. 



5 How does non-linear phase determine amplitude? 

In this section, we shall completely characterize a real- valued function p G iSa. We begin with 
introducing two one-sided filters that are related to the two-sided symmetric cascade filter Ha'- 

H+{t) := /7,(t)xM+(t), H-{t) := i7,(t)xR_ (t), (5.1) 

where t G M. The next lemma provides the Fourier transform pairs of and , respectively. 



Lemma 5.1 The Fourier transform of and are given, respectively, by 

1 1 1 - e-'« 



and 



27r 1 - ae-i€ i^ 



1 1 - e*^ 



(5.2) 



1 



-i^ 
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Proof: We start with the observation 

kez+ keN 
By the definition of the Fourier transform ( \1.3\\ . we obtain that 

1 — a 



V2tt 
1 - a 



"^■^ fcGN 



fcGN 

e dt 



V27r J[o,fc) 

where, the interchange of the order of the summation and the integral is justified by the absolute 
convergence of the series. Continue by noting XlfceN '^'^^^^''^^ ~ i~a&i have 

W(-)r(f) = ' 1 — - (1 - a)e-^ 



^ 1^(1 - ae-i«) 
1 1 l-e-'« 
1 - i^ 

Using the identity H~{-) = H^{—-), we obtain that 

1 1 1 - e^« 



/27r 1 - ae'« -i^ 



□ 



Lemma 5.2 The quasi- Bedrosian type identity 



H (pa(-)sinc(-)) (t) = i^p„(t)?{sinc(t) (5.3) 

1 — a 



zs trae /or t G M and 

1 pi* - 1 1 

(sinca(t) + iHsinCa(t)) . (5.4) 



1 — ae'* it 1 + a 

Proof: Let r(t) := yz^^^it^- Separating the real part from the imaginary part of r{t) produces 

1 — a sint 1 + a 1 — cos t 

rit) = h 1 

1 — 2acost + a^ t 1 — 2acost + a^ t 

1 1 

Pa(t)sinc(t) + i- pa(t)'Hsmc{t) , (5.5) 



1 + a 1 — a 
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where we have used the identity 

^ , . /X 1 — cost ,^ ^, 

'Hsinc(t) = . (5.6) 

On the other hand, observe that r{t) = y/2TT (H^)^ {~t), so the Fourier transform r{t) = 
\/27iH^{t) has zero negative spectrum. This imphes r is an analytic signal. Therefore the 
imaginary part of r{t) equals to the Hilbert transform of its real part, that is, 

n ( -^pa(-)sinc(-) ) (t) = --^Pa{t)'Hsmc{t). 

After rearranging the above equation we obtain equation (15. 3p . Now we turn to show equation 
(15. 4p . Applying equation (15. 3p to the imaginary part of equation (15. 5p we continue to have that 

r{t) = ^(p,(t)sinc(t) + iH(pa(-)sinc(-)))(t) 

= T^— (sinca(t) + i'HsinCa(t)) , 
1 + a 

which is equation (15.41) . □ 

Corollary 5.3 Let cosinca(t) := and define cosinCa (0) = lim^^o then for 

t G M, 

Hsmcait) = cosinca(t). (5.7) 

Moreover, the function cosine^ is odd, bounded, infinitely differentiable, cosinCa(t) < (yt^)^ 
for t G R, and cosine^ G L^(R). 

Proof: Equation (15. 7p is just a new form of equation (15. 3p . Indeed, the left-hand side of equation 
(15. Sp is just 'HsinCa(t). Utilizing equations (11.70 and (15. 9p then the right-hand side of equation 
(15. 3p is simplified to cosinca(t). The second statement follows by rewriting cosinCa as 

/ \ 1 + a / \ 1 — cost ,_ ^, 

cosmca(t) = pa[t) , (5.8) 

1 — a t 

and noting the function 

^ , . /X 1 — cost 

7{sinc(t) = ^ (5.9) 

is infinitely differentiable if we define 'Hsinc(O) = lim^^o^— and [ [ < ^-i^ for any 
t G M. □ 

An example of the graph of cosine^ with a = 0.5 is shown in figure 15^21 As a comparison, 
the graph of 'Hsinc(t) = i^^^^ is also shown in the figure. 

The following lemma appeared in |3]. However, a completely new proof by direct construc- 
tion is given here for both sufficiency and necessity. The construction approach casts a new 
insight to understand the relation of the spectrum p of a function p G 5a to the symmetric 
cascade filter Ha. 
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Figure 5.2: The graph of cosinCa, a = 0.5. 

Lemma 5.4 A real signal p ^ Sa if and only if there are two real sequences r = {rk : A; G Z} 
and s = {sk : A; G Z} m /^(Z) such that for t G M, 

p(t) = ^rfcsmca(t-2fc7r) + ^Sfc('HsmCa)(t-2A;7r). (5.10) 
feez feez 

Proof: We first show the necessity, that is, if p G iS^, then p is given by equation fl5.10l) . From 
Lemma IXBl we know that p(OxiR+(0 p(OxiR-(0 completely determined by p(OX(o,i)(0 
and p(OX(-i,o)(05 respectively. Expand p(OX(o,i)(0 p(OX(-i,o)(0 ii^to 1-periodic functions 
7i(0 := ^Mfce'^''''* and 72(0 := ^^fce'^'"'*, where the two series {uk : k e Z) and {vk : k e Z) 

are in P{Z). 

By Lemma [3.61 we thus obtain that for a.e. ^ G M, 

p(0 = 7i(0^:(0 + 72(0^-(0- 

The Hermitian property of p implies that for ^ G R \ {0}, 

72(o = (7i(-e)r. 

For convenience let 7 = 71. Consequently p has the representation 

for a.e. ^ G R. Let pA be the analytic signal associated with p, which may be defined by 

p^(0 = 2p(0xM+(0 
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or equivalently for t G M, 

PA{t)= p{t) + mp{t). 

Then we have p{t) = Re(pyi(t)). We need to investigate pA- The inverse Fourier transform of 
Pa yields 



PA{t) = ^ [ 27(0^a^(0e^^Me 



Applying the Lebesgue dominated convergence theorem to interchange the order of the integral 
and the sum, and appealing to equation fl5.2p we obtain that 



2 I gi(t-2fc^) _ 2 



/2^^ ''l-aei(*-2M Ut-2kn 



Denote Uk = u^;}^ + iu^^\ Recalling equation fl5.4l) we readily obtain that 

P{t)- 



E4'^sinc,(t - 2k'K) - (^sinc,(-))(t - 2kn) ] . 



'2ti{1 + a) 

Let Th = r-^,, — tUi' and Sh = — 7=7; — tU^u^ . Since the series (uk,k G Z) G /^(Z), we must 
have both series (r^, A; G Z) and (s^, G Z) in /^(Z). We have arrived at equation flS.lOp . 

We now turn to the proof of sufficiency. It suffices to check that a function p G L^(]R) having 
the form (15.101) satisfies the equation f l3.30p by Lemma 13.61 Applying the Fourier transform to 
both sides of equation flS.lUI) and noting that the Fourier transform of the generalized sampling 
function given by equation (12.170 we get that 

m = HaiO {M{0 - i sgn(e)G(0) 
with the two 1-periodic functions 



M(0 = ^(1 + a) 5^r,e-'2^-«, G(0 = ^(1 + a) ..e^^'^-?. 

fcez fcez 

Observing the identity Ha{- + n) = a^'^^Ha{-) and the 1-periodicity of M and G lead to that, for 

p{^+\n\) = H,(^+\n\){M{^+\n\)-iG{^+\n\)) 
= aH/7,(0(M(0-iG(0) = aHp(0. 

The case with negative ^ can be shown similarly. This completes the proof of the lemma. □ 
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Theorem 5.5 If p ^ Sa, then both p and Tip are continuous. 

Proof: By Proposition 13.51 Tip G Sa, thus it suffices to show p G 5a is continuous. Indeed, if 
p G Sa, then p has the representation f lS.lOp by Lemma [5.41 Since both sinca and "Hsinca are in 
L^(R) by Proposition 12.11 and Corollary 15. 3[ consequently by the Cauchy-Schwarz inequality, 
the two series on the right-hand side of equation f l5.10p converge uniformly, hence the limiting 
function p is continuous. 

□ 

Lemma 5.6 If f E Sa, then for k & Z, 

if, sinca(- - 2k7r)) = f{2kTi), and (/, 1-Lsmca{- - 2kTi)) = -l-Lf{2kTx). 
Proof: The 27r-periodicity of sin 9a yields 

(/,sinc,(- - 2kn)) = ^/(t)^^!^M_^dt = -?^(/sin^a)(2fc7r). 

Therefor by the assumption on /, we have when /c G Z, 

-H{f smea){2k'K) = /(2fc7r) cos0a(2fc7r) = /(2fc7r) 

because cos6'a(2A;7r) = cos^^a(O) = 1 by equation (I2.15p . 

The second equality follows by noting the fact that Tif G Sa due to Proposition 13.51 and the 
Hilbert transformer is an anti-self adjoint operator thus 

{f,nsmcai- - 2kTT) = -{nf, sincai- - 2kTx) = -nf{2k'n). 

□ 

Lemma 15. 4[ Lemma 15. 6[ and Theorem 12.21 immediately implies the following theorem. 
Theorem 5.7 Any real-valued function p G iS^ if and only if 

p{t) = 1^ \^p(2A;7r)sinCa(t - 2A;7r) - V 'Hp(2A;7r)?^sinCa(t - 2kTi). 

1 + a ^-^ 1 + a ^-^ 

Moreover, the sampling sequences (p(2/c7r), /c G Z) and {'Hp{2kiT), A; G Z) are in Z^(Z). 
From Theorem 15.71 and Theorem 12.21 we have the following corollary. 

Corollary 5.8 The system $ defined in f l2.2ip is a complete orthogonal system of the subspace 
Sa C L^(M). That is, the subspace equals to the closure of the span of the set $, or symbolically, 

Sa = span$. 
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Lastly, we state some facts for the case a = 0. Note when a = 0, sin 6a{t) becomes sint, 
cos9a(t) becomes cost, sinca becomes sine and "HsinCa becomes "Hsinc which is given at t G M 
by equation ( I5.9p . Therefore Theorem I4.H Theorem 15 . 71 and the well-known Shannon sampling 
theorem imply the following corollary. 

Corollary 5.9 The following statements are equivalent: 

1. The real signal p G is bandlimited such that suppp C [—1, 1]. 

2. n{p{-) cos(-))(t) = pit) sint, t G M. 

3. nipi-)e'-){t) = -ip(t)e", t G M. 

I Pit) = EkezPi^k7r)smc{t - 2kn) - Hp(2fc7r) ^^g^^'f t G M. 

5. p{t) = X]fcezP(^^)si^c(^ ~ ^^)' i £ I^- 
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